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Abstract 

We account for the different symmetries of the 2 1 ' 3 S' helium excited states in a quasiclassical 
description of the knockout mechanism augmented by a quantum shakeoff contribution. We are 
thus able to formulate the separate contribution of the knockout and shakeoff mechanisms for 
double photoionization for any excess energy from the 2 1 ' 3 S' states. Photoionization ratios and 
singly differential cross sections calculated for the 2 1,3 5 excited states of helium are found to be in 
very good agreement with recent theoretical results. 

PACS numbers: 32.80.Fb, 03.65.Sq 



1 



A two electron transition in an atom after the absorption of a single photon is only possible 
due to electron-electron correlations. The study of two electron transitions in helium, the 
simplest atomic target with two electrons, probes the role of electron-electron correlations 
in their purest form. As a result, there has been an extensive amount of experimental and 
theoretical studies of double photoionization from the ground state of the helium atom. 

In addition to the role of electron-electron correlations, the study of two electron escape 
by a single photon absorption from the singlet, 2 1 S, and triplet, 2 3 S, excited states of helium 
probes the role of symmetry in the two electron dynamics. This is particularly important for 
quasiclassical approximations where symmetries and quantum interferences are difficult to 
include. With the exception of a few earlier studies in the high energy limit 

HQ 

, it is only 

recently that the study of double photoionization from the 2 1,3 S' excited states of helium has 
attracted considerable theoretical interest. Experimental measurements of photon absorp- 
tion cross sections from the helium excited states still remain a challenge. In these recent 
stad.es sophisticated full. a„l at-haitio m ethods na m el y the R - m a«x m ethod □ aad 
the convergent close-coupling method |4] are used to compute integral double ionization 
cross sections, double to single ratios, and only very recently, single and triple differen- 
tial cross sections as well Q|. However, due to the numerical nature of these sophisticated 
ab-initio calculations approximate methods are also essential in uncovering the underlying 
mechanism of the two electron dynamics in the double escape process. Such approximate 
methods have been successfully developed for the double photoionization from the helium 
ground state . The current work focuses on an approximate method to describe double 
photoionization from the helium excited states. 

In the double photoionization process the redistribution of the energy, following the pho- 
ton absorption, is often discussed in terms of two mechanisms, that is, knockout and shakeoff 
dominant in the low and high energy limit, respectively. In the knockout mechanism one 
electron, the primary, absorbs the photon and undergoes a hard collision with the secondary 
electron, thus, knocking it out of the atom. The knockout mechanism dominates at low 
energies where the interaction time of the two electrons is large and can be described clas- 
S iea„ y QQ. For Knockout, elects eolations are essential to the anal eontinuu, state. 
On the other hand, in the shakeoff mechanism the primary electron absorbs the photon 
and leaves the atom very fast without undergoing a direct interaction with the secondary 
electron. The secondary electron feels the sudden change in the atomic field and relaxes in 
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one of the bound or continuum eigenstates of the remaining ion. The shakeoff mechanism 
is quantum mechanical in nature and prevails at high energies where the interaction time 
of the two electrons is short. For shakeoff, electron correlations are important in the initial 
state before the photon is absorbed. 

In the following, we formulate the explicit but separate contribution of both mechanisms 
for double photoionization from the 2 1 ' 3 S' helium excited states for any excess energy by 
extending the model developed to separate the contribution of both mechanisms for the 
helium ground state in ref.p]. The separation is achieved by a quasiclassical formulation 
of the knockout process. By construction it is free from any shakeoff contribution which is 
purely quantum mechanical. This separation not only facilitates the calculation of double 
photoionization but offers considerable insight into the process concerning for example the 
similarity with electron impact ionization of He + . Compared to the ground state, formulat- 
ing the separate contribution of both mechanisms from the helium excited states is a much 
harder problem since the different symmetry of the 2 1 ' 3 S' states has to be accounted for in 
the model. 

We express the two steps of the double photoionization process from the helium excited 
states, absorption of the photon and afterwards redistribution of the energy, as 

<4 + = a abs P+ + , (1) 

where X stands for either knockout or shakeoff and P£ + is the double electron ionization 
probability. In what follows we discuss how to obtain P^o an d then P$q for the 2 1 ' 3 S' 
helium states. 

-Knockout mechanism- After the photon absorption by the primary electron, we describe 
the subsequent evolution of the two electrons using the classical-trajectory Monte Carlo 
(CTMC) phase space method. CTMC has been successfully used to describe charged particle 
impact processes We model the initial phase space of the trajectories 



quantum mechanically using a Wigner distribution [lj] . To do so, we first find the initial two- 
electron wavefunction. Since electron correlations are not important in the initial state of the 
knockout mechanism, we choose our initial state as the independent electron wavefunction 

^o(ri,r 2 ) = 0f-(r 1 )0ff(r 2 ), (2) 

where ^(n) = (Z^/vr) 1 /^-^ 1 and (r 2 ) = (2| b /(8tt)) 1 / 2 (1 - Z 2s r 2 /2)e- z '^/ 2 are 
hydrogenic Is and 2s orbitals with effective charges Z ls and Z 2s , respectively. ty (ri,r 2 ) 



given by Eq.(J2J) is not a symmetrized wavefunction, thus, we do not yet account for the 
different symmetries of the triplet 2 3 S and the singlet 2 l S helium excited states. To find 
the effective charges Z la and Z 2s we use the two-electron ionization energies I 2 i + s / 2 z s ~ 
2. 146/2. f 75 for the 2 1 S/2 3 S states given in ref.Q, and the ionization energies of the Is and 
2s orbitals as follows. The energy needed to remove the Is electron after the 2s electron has 
been ionized is 2 a.u. and thus from Ii s = Zf s /2 = 2 we find Z\ s = 2 for both the triplet 
and the singlet states. From I 2 i%/ 2 -i S ~ -^ ls = ^ 2s = Z 2s /% we find Z 2s ~ 1.081/1.184 for the 
2 1 S/2 3 S states, respectively. Atomic units are used throughout this paper unless otherwise 
indicated. 

While, for the Is 2 helium ground state the photon is necessarily absorbed by a Is electron, 
for the ls2s configuration of the 2 1 ' 3 5' helium states the photon can be absorbed by a Is or a 
2s electron. For the values of the effective charges given above, one can show that the cross 
section is much larger for photon absorption from a Is electron rather than a 2s electron 
using the independent electron picture [la]. Therefore, we take the primary electron to be 
the one on the Is orbital. With the PEAK approximation, exact in the high photon energy 
limit Q], we assume that the photon absorption happens directly at the nucleus. This 
significantly reduces the initial phase space to be sampled by the CTMC method. 

The initial phase space distribution, p(T), is the two electron density immediately after 
the photon absorption given by 

p(T) = N5(r 1 )p 2 (r 2 ,p 2 ), (3) 

where N is a normalization constant. The initial distribution of the primary electron is <5(Vi) 
(PEAK approximation) while that of the secondary electron is given by 

P2(r 2 ,p 2 ) = W^r 2 ,p 2 )5(et l -e B ), (4) 

where W^(r 2 ,p 2 ) is the Wigner distribution function of the two electron wavefunction with 
the primary electron at the nucleus, r\ = 0, 

1>(r 2 ) = * (n = 0,r 2 )((* (ri = 0,r 2 )|* (ri = 0,r 2 )))^ 2 . (5) 

In Eq.(J3J), we take the energy of the secondary electron immediately after photon absorption, 
e| n , to be fixed on the 2s energy shell e# = —Z 2s /8. From Eq.(J2J), it follows that e l 2 = 
pl/2 — Z 2s /r 2 . The excess energy available to the two electron system after photon absorption 



is determined by the photon energy to be 



E = u-I+% E = u-I+ + S (6) 

for the 2 1 S and 2 3 S states, respectively. Due to the PEAK approximation the primary 
electron can have any energy necessary so that together with the initial energy e?> n of the 
secondary electron it adds up to the excess energy E in Eq.flBJ). After modelling the initial 
phase space distribution, we propagate the electron trajectories using the classical equations 
of motion (CTMC). Regularized coordinates are used for the propagation of the electron 
trajectories to avoid problems with trajectories starting at the nucleus (t*i = 0). Doubly 
ionized are those trajectories that end with the asymptotic energies of both electrons being 
positive. To evaluate the double electron escape probability each trajectory is weighted by 
the initial phase space Wigner distribution. 

So far we have treated the two electrons as distinguishable particles, that is, we dis- 
tinguish between the primary and the secondary electron. To account for the singlet and 
triplet symmetries we have to symmetrize the probability amplitude (differential probability) 
dP^o/dc with respect to the two identical particles, 

dP& _ 1 UdP^E) ± ^++(£-e,£) ^ (7) 

In Eq.(J7J), dPftQ^e, E)/de is the probability for both electrons to escape when the primary 
electron is ejected with energy e, where < e < E, and the secondary electron is ejected 
with energy E — e. To evaluate dP^Q^e, E) / de we divide the energy interval [0, e] into N 
equally sized bins and find the doubly ionized trajectories which fall into the bins. In our 
calculations we take N = 21 for excess energies up to 80 eV and N = 27 for higher excess 
energies. The double ionization probability P^o ls obtained by integrating over all possible 
energies that an electron can be ejected with, that is 

Note that for the case of the helium ground state the double ionization probability P^o ls 
worked out without using the differential probabilities f|. 

-Shakeoff mechanism- Assuming that the primary electron is suddenly removed from the 
atom, Aberg found that the probability for the shaken (secondary) electron to relax on 
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a hydrogenic eigenstate of the remaining ion for any excess energy is 

p: = mw)\ 2 i(rw), (9) 

where ip u {r 2 ) = J <i 3 r 1 z/*(r 1 )\I/o( r *i) r 2 ) and v{v\) is the primary electron wavefunction after 
it has left the atom. The primary electron is in an s state before the photon absorption and 
in a p state afterwards. x I / o( r 'i ; r 2) is the initial state wavefuntion of the 2 1 ' 3 S' helium excited 
states and (fi a is a hydrogenic eigenstate of the bare nucleus [Z = 2) that is either a bound 
(a = n) or a continuum state (a = e). Aberg has further shown that when the primary 
electron leaves the atom with very high energy (^(t*i) = (2n)^ 3 ^ 2 e^ l ^ iri ) Eq.(JH} takes the 
simplified form 

p \M%(ri = 0,r 2 ))| 2 

a (*o(r 1 = 0,r a )|* (ri = 0,r 2 )>- 1 J 

Eq. (|10j) reveals the quantum mechanical nature of the shakeoff process since it is expressed 

as an overlap of the initial bound state wavefunction and the final continuum state wave- 
function. 

Although Eq. (jlUj) was derived in the high energy limit we assume that the primary 
electron absorbs the photon on the nucleus for all excess energies, that is, we adopt the 
PEAK approximation as in the knockout case. To find the double escape probability PgQ 
we then integrate over all possible energies of the shaken electron in the continuum 

Pto(E)= [ E pje. (11) 
J o 

We further simplify the evaluation of the shakeoff probability in Eqs. (jlU)) . (jllj) by taking 
the initial state to be the symmetrized wavefunction 

Mri,r 2 ) = Nx ($ {Tx)4$ {T2)±4$ {r2)4$ {rx)\ (12) 



for the singlet and triplet states, respectively, with N\ a normalization constant. The initial 
state correlations are accounted for only through the effective charges. We next assign the 
same set of effective charges Z\ Q and Z\ Q for both the triplet and the singlet states as 
follows. The asymptotic ratio (high energy limit) of double to single ionization is found very 
accurately in ref.|2:] to be = 0.009033/0.003118 for the singlet /triplet states where in 
our model Rqq is given by 

Rao = P$${E - oo)/(l - P++(E - oo)), (13) 



and 

P${E oc) = / P € de = l-J2 Pn- (14) 
Jo 

Using Eqs.([13|). ()14j) and the symmetrized wavefunctions in Eq. (jl2}) we find the sets of 
charges that match both asymptotic ratios for the singlet and the triplet states. We 
then select that set of charges for which the shakeoff double ionization probability as a 
function of the excess energy, obtained using the simple wavefunctions given in Eq. (fT2l, i s 
closest to the one obtained using the fully correlated Hylleraas wavefunctions given in ref . 19( . 
The set of charges thus found is Z x so « 1.757 and Z\ Q ~ 1.728. The reason we do not use 
the Hylleraas wavefunctions given in ref.[l9] is that they do not reproduce the accurate 



asymptotic ratios obtained in ref. [2| using highly accurate Pekeris-type wavefunctions. We 
emphasise though that one does not need to use the approximate wavefunctions in Eq. (jl2j) 
to compute the double ionization probability; highly accurate wavefunctions that reproduce 
the correct asymptotic ratios could be used instead. 

For the shakeoff probability, P e in Eq. (jl(Jj) already gives the differential double ionization 
probability. Despite the symmetrization in Eq.([12p we have lost the indistinguishability of 
the electrons by identifying one electron as the primary one which absorbs the photon. Thus 
we need to symmetrize again in the final state with respect to the equal energy sharing point 
e = E-e = E/2 [6]. That is, 

dP& l(dP e dPB-A n ,, 

-Photoionization ratios-. According to Eq.([T} a ++ = ct o6s (P^q + Pso) an d o" + = cr abs — 
a ++ . Thus, the double to single ionization ratio is given by 

= P++/(1_P++), (16) 

where P ++ = P^q + P^o ■ 111 figure (JTJ) we compare the double to single ratio for the 2 1 ' 3 S' 
helium excited states with the results obtained by Kheifets et al jj| using the convergent 
close-coupling method and show that there is a very good agreement. The agreement is 
better for the 2 3 5' state. We find that the deviation occurs, particularly for the 2 3 S state, at 
photon energies where the contribution of knockout and shakeoff mechanisms is comparable. 
At these energies any interference effect between the knockout and shakeoff mechanism would 
have its largest effect. So, it may be that the deviations we see are due to that interference 
effect that we do not account for in our calculation, since we add the knockout and shakeoff 



contributions incoherently. For the 2 1 S' state a maximum of ~ 2.84 % is reached at 14 eV 
above the double ionization threshold of the 2 1 S state. For the 2 3 S state a maximum of ~ 
0.69 % is reached at 60 eV above the ionization threshold of the 2 3 S state. In figure 
we see that at high energies the knockout contribution goes to zero as expected and the 
shakeoff contribution dominates and reaches the asymptotic limit of 0.009033/0.003118 for 
the singlet /triplet states. Stronger correlation effects for the singlet symmetry (ri = r<i is not 
forbidden as is the case for the triplet) result in a much higher double to single ionization 

ft 

ratio compared to the triplet case. From figure (JH), for the 2 1,3 S' states, and ref. p, for 
the helium ground state 1 S, we see that as we go from 1 S — > 2 1 S — > 2 3 S the shakeoff 
mechanism overtakes the knockout mechanism at smaller energies. The reason is that the 
electron-electron correlation becomes smaller, thus diminishing the knockout contribution 
and favouring the shakeoff mechanism at even smaller energies. 

-Single differentials- To compute the single differential probabilities for the helium excited 
states we use 



where for a a b s we use the results for the total-photoionization cross section for the 2 5 and 
the 2 3 5* states given in ref.Q. In figure (j2J) we compare our results for da ++ /de with the 
results obtained very recently by Colgan et al \^ using the convergent close-coupling method 
for four values of the excess energy. We see that our results for the single differential cross 
section as a function of the ejected electron energy normalized by the excess energy are 
smaller for the 2 1 S' state while there is an excellent agreement for the 2 3 S' state. Again as 
we go from 1 S — > 2 1 S' — > 2 3 S' the single differential cross sections are more U-shaped for the 
same excess energy. The reason is again that the electron-electron correlation decreases thus 
favouring the ejection of one fast and one slow electron, that is favouring unequal energy 
sharing. 

In figure (jSJ) we show the separate contribution of knockout and shakeoff to the single 
differential probabilities for the 2 1,3 S' helium states for excess energies 10 eV, 40 eV and 160 
eV. For the singlet state the knockout contribution dominates at small excess energies, 10 eV 
and 40 eV, while as the excess energy is increased to 160 eV the shakeoff contribution begins 



de de de 

In addition we compute the single differential cross sections using 




(17) 




(18) 
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to dominate regions of unequal energy sharing. For the triplet case the shakeoff mechanism 
is already significant at small excess energies. Note, that the knockout contribution for the 
triplet case is zero at the equal energy sharing point, e = E— e, because of the symmetrization 
with respect to the two identical electrons, see Eq.(J7j). 

In conclusion we have shown that the double ionization from the 2 1,3 S' states can be 
accurately described by a separate formulation and calculation of the knockout and shakeoff 
mechanism at any excess energy. In comparison to the helium ground state jo| this is a 
harder problem because we have to account for the different symmetries of the singlet and 
triplet states. The success of this simple model to describe double ionization from the helium 
ground state as well as the helium excited states is proof for its validity. In the future, we 
plan to use this simple model to describe triple photoionization cross sections. 

The authors wish to thank T. Pattard for helpful discussions. 
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List of Figures 



• Figure 1. Double to single ionization ratio as a function of the photon energy. 

n 

Dots/open circles indicate the results of Kheifets et al [4] in the velocity/acceleration 
gauge. For the triplet state Kheifets results in both gauges are indistinguishable to the 
scale of the figure. Our results are indicated by a solid line for the total, by a dashed 
line for the knockout and by a dashed-dot line for the shakeoff double to single ratio. 

• Figure 2. Absolute single differential cross sections as a function of the electron ejected 
energy scaled by the excess energy. The dashed lines are the results by Colgan et al. 
Our results are indicated by solid lines. 

• Figure 3. Absolute single differential probabilities as a function of the ejected electron 
energy. The knockout contribution is indicated by dashed lines while the shakeoff by 
solid lines. The top panel is for the 2 l S state while the bottom is for the 2 3 S state. 
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